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Recent experiments on particle diffusion in bacterial baths indicate the formation of cor- 
related structures in the form of bacterial swirls. Here we predict that such a structural 
ordering would give rise to the new effect of propagating shear waves in a bacterial bath 
at length scales of the order of a swirl, which corresponds to time scales of the order 
of the lifetime of a swirl. Our prediction can be tested against future experiments in 
bacterial baths. 
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1. Introduction 

A simple liquid cannot sustain shear stress under a static condition. However, at 
molecular length and time scales (i.e. time scales of the order of the relaxation time 
of the liquid), it is known that a liquid behaves much like a solid with a viscoelastic 
character pp. One can understand such a phenomenon qualitatively as follows. On 
molecular length scales a dense liquid develops solid-like structures due to slowing 
down of diffusional structural relaxation and consequently it can support a shear 
wave. Molecular dynamics simulations of hard sphere liquids and experiments on 
liquid Argon support this picture |2I3| . On the theoretical end, such a behavior can 
be understood in terms of viscoelastic models 

In recent years there have been experiments probing 'Brownian Motion' of 
polysterene spheres in bacterial baths ;5;. These experiments point to a crossover 
from superdiffusive to diffusive behavior of the mean-square displacement of the 
polysterene beads suspended in a bacterial bath. This effect is similar to ballistic 
to diffusive transition in colloids 0E1IZ] where the borderline for such a transition 
is given by the viscous damping time to = ^ (m& being the mass of the Brownian 
particle and 7 the friction coefficient), which for a polysterene bead of diameter 
lO/im is about 10 _5 sec. However, in the case of a bacterial bath one notices that 
such a transition takes place at a time of about r c w lsec. This characteristic time 
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scale r c increases with the number density n of bacterial cells in the bath. The 
length scale l c — y/[Ar 2 (t c )], which is the root mean square displacement at t = t c 
also increases with n. These observations point to the formation of correlated struc- 
tures in the form of bacterial swirls which have been observed in these systems. In 
fact, one finds that l c is indeed of the order of the spatial scales of these swirls 
(10 — 20/im). To sum up, one of the key observations coming out of these exper- 
imental studies is the emergence of a time scale and a corresponding length scale 
in the Brownian motion of a polysterene bead suspended in a bacterial bath. This 
time scale r c increases with the number density n of the bacterial bath, indicating 
the appearance of solid-like ordering over length scales around l c = ^/[Ar 2 (r c )]. In 
a dense liquid such ordering takes place on molecular length scales. In the present 
experiments the bacterial bath can be viewed as a scaled up liquid with a typical 
relaxation time scale w lsec in contrast to a colloidal liquid characterized by a time 
scale w 10~ 5 sec. 



2. Viscoelasticity in Bacterial Bath 

We now describe in some more quantitative detail the crossover from super-diffusive 
to diffusive behavior in a bacterial bath. Consider the dynamics of a bead of mass 
mi, suspended in a bacterial bath j^j in terms of a phenomenological Langevin-type 
equation: 

m b — = --fv + f(t) (1) 

where 7 is the friction coefficient of the fluid. The collisional force f(t) experi- 
enced by the bead due to "kicks" from the swirls of bacteria has a zero mean and is 
exponentially correlated in time with a time scale r c corresponding to the lifetime 
of a swirl. Thus 

< /(<) >= 0; < /(0) • f(t) >=< / 2 (0) > ex P (-t/ Tc ). (2) 

where < / 2 (0) >= Dj 2 /t c with D = ksT/j the Einstein diffusion coefficient. Such 
an analysis leads to the following behavior for the mean square displacement of 
the bead. At short times (t « t c ) the motion of the bead is ballisticlike with 
< Ar 2 (t) >~ t 2 whereas at long times (t » t c ) it is diffusive < Ar 2 (t) >~ t 0. 

More recent experiments on correlated bacterial motion also indicate the 
emergence of a long time scale t c increasing with the number density n of bacteria. 
One interesting observation that comes out of the study jH] is an enhancement of the 
viscosity 77 of the bacterial bath. This enhanced viscosity is proportional to the large 
time scale r c emerging in the problem. In this paper we show that this crossover 
time scale t c is connected to a viscous to an elastic transition in a bacterial bath. 
1 Let us consider the following model which interpolates between a purely viscous 
response and a purely elastic one for a liquid of viscosity 77 and high frequency 
rigidity modulus G. 

,1 1 3, dv x dv z . . 

y~ + T< Tj a xz = IT + IT ■ ( 3 ) 
77 G at oz ox 

J We emphasize that the main aspect of a bacterial bath that we are interested in is the presence of 
correlated structures in the form of swirls. Such structures give rise to a large spatial scale in the 
problem. This in turn leads to a large crossover time scale t c which is of the order of the lifetime 
of a swirl. 
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Here a xz is the xz component of the shearing stress tensor corresponding to the 
strain rate + . On taking the Fourier-Laplace transform of Eq. one notices 
that the viscoelastic approximation corresponds to replacing the zero frequency 
shear viscosity rj by a frequency dependent shear viscosity: 

*?H = ( ■ ~T g \ ' ( 4 ) 
Defining r c = ^ to be the relaxation time, we can rewrite Eq. Q as 

m = , . G M i V (5) 

Therefore on long time scales for which lut c <C 1 one expects the flow to be viscosity 
(r) = Gt c ) dominated. As one starts probing shorter time scales, one expects a tran- 
sition from a viscous response to an elastic one to take place for lut c > 1. Since for 
a bacterial bath r c w lsec we arrive at the following important quantitative point. 
Unlike a colloidal bath where the onset of propagating shear waves takes place on 
time scales w 10~ 5 sec, in a bacterial bath we predict that this would happen on time 
scales « lsec. Thus, if a bacterial bath is perturbed at a frequency of ks 1Hz one 
would notice a transition from a viscous relaxing behavior to an elastic propagating 
one. Piezoelectric bimorphs can be used to produce and detect controlled shear 
motions at a frequency ~ 1Hz jl(J| in a bacterial bath. In contrast, in order to 
see the propagation of shear waves in a colloidal liquid one needs to perturb it 
at a high frequency » 10 5 iJz. A bacterial bath thus provides us with a "scaled 
up" non-equilibrium statistical mechanical system. This is the key qualitatively 
and quantitatively interesting observation made in this paper. We will analyze the 
emergence of shear waves in some more detail below. 

An analysis based on memory function formalism can be used to obtain a sta- 
tistical mechanical description of the transverse current fluctuations in a bacterial 
bath. In this formalism the spectrum CT{k,cu) of transverse current fluctuations ( 
i.e. the Fourier-Laplace transform of the transverse current autocorrelation function 
C T (k,t) =< ji(t)j^U0) > ) takes the form: 

C (k w) = = f 6 ) 

Hw + £ M r (fc,w)]' 

where Mr(fc,w) is the pertinent memory function which converges to the hydro- 
dynamic viscosity r] in the long wavelength and long time limit. Modelling the 
memory function in the time domain as Mx(k, t) = Mx(k, t — 0)e~ t / Tc ( fe ) we arrive 
at the following expression for the real part CT{k,uj) of the spectrum CT(k,uj) of 
transverse current fluctuations: 



uj 2 ujl t T c (k) 

7r[^+T2(fc)K-^)2 



where r c (k) is the wave vector dependent lifetime of correlated structures in the 
bacterial bath; ujq — Cr{k,t — 0) and uj\ t — £-SL^ with p the mass density of a 
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bacterial swirl. We can now establish the criterion for the existence of propagating 
transverse modes. From Eq.Q it is clear that the condition for Cx(k,uj) to have 
a peak at a nonzero frequency, i.e. the threshold condition for shear waves to 
propagate at a finite speed is ujf t T^(k) = i. The corresponding value for the 
critical wavevector is k c = j 1 - — -^-\f^ which corresponds to a length scale of 
the order of the spatial scale l c of a bacterial swirl. Thus, as one starts sweeping 
over wavevectors, the form of the spectrum Cy(fc, ui) changes from a profile peaked 
around lot c — for kl c << 1 to one peaked around a finite nonzero frequency 
ujt c « 1 for kl c « 1 (See Fig. 1). 
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Fig 1. Schematic representation of spectrum of transverse current autocorrelation function 
Cr(fc,^) as a function of ojt c . The figure shows transition from a profile centered around lut c = 
(dashed curve) for kl c << 1 to a profile with a side peak at uit c = 1 for kl c = 1 illustrating the 
phenomenon of onset of shear waves. 



Using r c pa 2sec, l c w lOfim and an appropriate value of p based on experimental 
observations |S] one can thus calculate the high frequency rigidity modulus G for a 
bacterial bath. 



3. Discussion 

There have been some interesting studies of Brownian motion of spherical probes in 
wormlike micelles where they study the role of polymer solutions in determining 
diffusion behavior of probes. The focus of the work presented here is quite different 
from that. We focus on a super-diffusive to diffusive behavior of the type that is 
observed, for instance, in a simple liquid where a ballistic to diffusive transition 
takes place at a typical relaxation time scale This leads to an interesting 

viscoelastic effect. To summarize, we predict an onset of shear waves in a bacterial 
bath at rj 1Hz as a result of structural ordering on larger length scales compared 
to that in a dense viscous colloidal liquid. This is a consequence of a large spatial 
scale arising from the swirl size of the bacteria forming the bath. This enables us 
to view a bacterial bath as a 'scaled up' dense viscous liquid. Furthermore, our 
work points to an independent test of structural ordering in a bacterial bath. We 
expect this work to generate interest in testing the predictions through computer 
simulations and experiments in bacterial baths. 
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